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I. INTRODUCTION 

Let A4 n be the set of complex n x n matrices and let \\T\\ denote the spectral norm of 
T € Ad n . We denote by a = cr(T) the spectrum of T and by uit its minimal polynomial. 
We denote by rriT the degree of my. In this article we are interested in estimates of the 
type 


\m,T)\\ <<&(|mr|,<T,C) (1) 

where R (£, T) = (£ — T) _1 denotes the resolvent of T at point £ G C — cr and $ is a 
function of \rriT\, cr and £. We present a general method that yields optimal estimates 
of the above type and allows us to pinpoint the “worst” matrices for this problem. As 
it turns out for any £ 0 cr among these matrices are so-called analytic Toeplitz matrices. 
Generally, Toeplitz matrices are characterized by the existence of a sequence of complex 
numbers a = (a/ c )^Z“~ 1 |_ 1 with 


T = 


ao 

«i 


\ On—1 


0-1 


ai 


O—n+l \ 


0-1 

ao / 


We call T a an analytic Toeplitz matrix if ak = 0 for all k < 0 and we will denote by 
T n CZ Ad yi the set of Toeplitz matrices and by '/,, CZ ‘ 1 the analytic T'oephtz matiices. 

We will think of an estimate of the above type as an extremal problem. That is, 
under certain constraints on the set of admissible matrices T and fixed £, we are looking to 
maximize || R(£, T)|| over this set. To guarantee a finite bound we at least need a constraint 
on the norm and the spectrum of possible T. In our discussion we will generally impose that 
imi < 1 and call such T a contraction. Note that this condition can always be achieved 
by normalization T / \\T \\. C n C A4 n denotes the set of contractions. To ensure that the 
resolvent is finite, ( must be separated from cr. Depending on the situation it will be 
convenient to measure this separation in Euclidean distance 


d(z, w) := \z — w\ , z, w £ C 


or pseudo-hyperbolic distance 


p(z, w ) : = 


z — w 


1 — zw 


, z, w £ C. 


The case ( = 0 in (1) corresponds to the study of the condition number CN(T ) = 
\\T\\ ■ ||T -1 ||. It was shown by Kronecker in the 19th century that 


1 


sup {|T 1 1 : T eC n , p{ 0, a) = r] = — 


( 2 ) 
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The upper bound can be seen using polar decomposition [NN1, p. 137]. To prove the 
reverse inequality, one can use A. Horn’s theorem (see [AH] and also [MO, Chap. 9, Sect. 
E]) that gives a converse to the famous Weyl inequalities [NN1, p. 137]. 

The case |£| = 1 was studied by Davies and Simon and applied to characterize the 
localization of zeros of random orthogonal polynomials on the unit circle [DS]. The corre¬ 
sponding estimate for |£| = 1 is 

sup{d(C, cr) ||i?(C, T)|| : T G C n } = cot(-^-). (3) 

Note that |£| = 1 implies that we have p(£, a) = 1 for any a such that we could add this 
condition to the optimization problem. In this paper we provide an approach, which 

1) allows us to treat the problems (2) and (3) simultaneously and generalize them to 
any |C| < 1, 

2) strengthens the result (3) for |£| = 1, by considering the sup over T G C n with fixed 
minimal polynomial and 

3) provides explicit analytic Toeplitz matrices that achieve the estimates (1), (2) and 
(3). 

Already for ( = 0 it is a non-trivial question which matrices achieve the estimate 
(2). A theoretic classification of such matrices was given in [NN1], However, it seems 
rather difficult to compute explicit matrix representations from this classification. It was 
suspected in [RZ] that the estimate (2) is not reached by any Toeplitz matrix; this claim 
is refuted here. 


II. RESOLVENT BOUNDS 


Before we state our main results we need to introduce some further notation. We denote 
by ID = {z G C : \z\ < 1} the open unit disc in the complex plane, B is its closure and <9B 
its boundary. For any finite set ffCDwe define 


s(C, cr) := max 


1 ~ |A| 2 
|1 —AC| 


A G cr, p( C, A) = p(C, cr) 


The function s is related to the Stolz angle [JC, p. 23] between (gD-ct and A G cr. For 
T G C n the function Q H > s(£, a(T )) is bounded on B by 1 + p(T ) < 2, where 

p(T) := nrax{|A| : A G cr(T)} 


is the spectral radius of T. Of particular importance to our discussion will be the analytic 
Toeplitz matrix X r p given entry-wise by 




0 

j.n—1 

p r n-{i-j+ 1 ) 


if i < j 
if i=j 
if i > j 


( 4 ) 


with r G [0,1] and /3 G [0,2], The spectral norm of X r> p is computed below in Proposi¬ 
tion II.6. 
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A. Optimal upper bounds 

Theorem II.1. Let. T G C n with minimal polynomial m and spectrum, a. For any £ G D— a, 
it holds for the resolvent of T that 

< 5 > 

where X r ^ is the \m\ x \m\ analytic Toeplitz matrix defined in (4), r = p(f, a) and (3 = 
s(C, a). Moreover, for any £ € [—1, 1] and A G (—1, 1), A ^ ( the analytic n X n Toeplitz 
matrix 


rj-\* _ 


A 0 

1 - A 2 A 

—A(1 — A 2 ) 1 — A 2 


o\ 


A 


: •• •• •• 0 

V (-A) n_2 (l - A 2 ) ... -A(l - A 2 ) 1 - A 2 A/ 

is a contraction of minimal polynomial m = (z — X) n that achieves (5). 

Remarks: 


1. Setting ( = 0 in Theorem II. 1 we recover Kronecker’s result for condition numbers. 
Moreover, an analytic Toeplitz matrix is extremal for the Kronecker maximization 
problem (2). This answers the questions raised in [RZ]. 

2. Considering Theorem II. 1 for £ G (90, and bounding 

||^1,s(C,<t(T))|| < l^i, l+p{T) | < ||^i,2|| , 

we recover the resolvent estimate of [DS] as = cot(^) (cf. Proposition II.6 

below). However, if e.g. p{T) is known we can improve on [DS]. 

3. Equality is achieved by a so-called model matrix corresponding to a fully degenerate 
spectrum. This matrix arises as a representation of the compressed backward shift 
operator of H 2 with respect to the Malmquist- Walsh basis , see Section III below for 
details. The model matrix is a natural generalization of a simple Jordan block of 
size n to the situation when the spectrum is a = {a} C B. 

The next proposition treats the case cr(T) C <9D. 

Proposition II.2. Let T G C n with spectrum, a C c®. Then for any ( G C — cr 

1 

d{ C, o')' 


me, r)|| < 


( 6 ) 
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It follows that if T is a contraction such that a(T) C c® then it satisfies a Linear 
Resolvent Growth with constant 1, see [BN]. Clearly, any diagonal matrix in C n achieves 
the inequality (6). 

Theorem II. 1 has the disadvantage that the upper bound depends on the slightly un¬ 
yielding quantity ||X r fl||. Taking together Theorem II.4 with Proposition II.6 (both below) 
leads directly to the following. 

Theorem II.3. Let T G C n with minim.al polynomial rn and spectrum cr. For any ( G D — a 
let p(C,cr) = r denote the pseudo-hyperbolic distance of Q to a. Then it holds for the 
resolvent of T that 

- (i(l,dC)rH(l- r |C|)' 

The estimate is sharp in the limit of large \m\. 


B. Constrained maximization of the resolvent 

For any given ( G D and r G [0,1] we study the quantity 

K((,r) := sup {d(l, a()\\R((,T)\\ : T G C n , p{(,a) = r} , (7) 

where o = {A : A G cr}. 1Z also depends on n but we keep the notation simple and do not 
write this explicitly. 

Theorem II.4. For any (gD and r G (0, 1) we have that 

^(C, r ) = —^ ll^r,/3maxll ) 

where /3 max = ■ The estimate is achieved by an analytic Toeplitz matrix and in the 

limit of large matrices we find 

rc(C, r )^^ 1 _ r | C | as n ^°°- 

Recall that Kronecker’s result (2) is 

K(0, r) = L 

Davies and Simon (3) have shown that 

nc, 1) = cot(^), C€9B-{1}, 

which can be seen in the above theorem by continuously moving f to <90. Note, however, 
the fundamental difference in the large n behavior of the resolvent for ( G D and ( G c®. 
This is a consequence of the fact that for £ inside the unit disk ||X ri/ 3 || is bounded from 
above, see Proposition II.6. 
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Finally for £ € 3D and a polynomial P a with zero set a so that p{<f, a) = 1 we can 
compute the constant 

/C(C, Pa) = supMC, <t) ||i?(C, T) II : re C n , m T = Pa} , (8) 

which is bounded by TZ((, r). 

Proposition II.5. Let P a be a polynomial with zero set a = a\ Ua 2 , ay C D and a 2 C 3D 
so that P a =P ai P (J2 . For any ( € 3D — a, we have 

/C(C, Pa) < max(l, j{X 1;S(C;(Tl) ||) < ||Xi j2 || = cot( p , ), (9) 

where X 1 , s ((,a)> -Xi, 2 € ,M|p |(C). Moreover, fixing f e 3D— a, ni, ?r 2 > 1 and 0 < p\ < 1 
and taking the supremum of /C(£, P CT ) ouer all sets a = ay U ay with n\ = \P ai \ and 
n 2 = \P<n\ such that maxAgo-j |A| < pi we get: 

SU P /C(C, Pa) = |^l,l+pi| (10) 

where X\,\ +pi e M ni (C). In particular, the supremum of the above quantity taken over 
0 < pi < 1 is equal to 


C. Computing ||X ri/ g|| 

In this subsection we study the quantity ||X rj( g|| for r e (0,1] and fi e [0, 2]. We provide 
an implicit formula for ||X r ^|. 

Proposition II.6. Let X Tj p with r G (0,1] and ft € [0,2] be the n X n matrix introduced 
in (4). The following assertions hold. 

i) If (3 = 1 — r 2 then |X rjj g|| = 1. 

ii) At fixed n, ||X r a\ grows m.onotonically with r and fi. 

Hi) At fixed r and 1 — r 2 < (3, |X rjj g|| grows monotonically with n. 

iv) Let U n denote the n-th degree Chebychev polynomial of second kind, i.e. U n (cos(6)) = 
blll sh"(e) )6,) • VP 7 ^ 1 — r 2 then |X r)j g|| = |A*|, where X* is the largest in magnitude 
solution of the polynomial equation 

p n+1 rU n + M n (A 2 - r 2n ~ 2 ((3 - l) 2 )t/ n _i (^j = 0 
with p. = A 2 + r 2n ~ 2 ((3 — 1) and 7 = —A 2 (r + 1/r) + r 2n ~ 2 (r + l(/3 — l) 2 ). 


( 11 ) 
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v) If r £ (0,1) and 1 — r 2 < f3 we have the limit 

| X r p | —> -j) n — * °°- 

vi) For any n, r £ (0,1), /3 max £ (0, 2] and 1 — ?’ 2 < f3 < /3 max we have 

\*r,li\ < 

The equations in iv) (see (11)) are discussed in Lemma II.7 below. For instance we can 
recover known values of ||Xi a\ from Proposition II.6, iv). In this case 

where 0* is the unique solution of (11) in [ 2 ”~ ] tt, tt) and it follows, see [OS], that 

|Vi,„|| = l. |X., 1 1 = 2sill( 1 _ s _ ) , IX 1 . 2 II =cot(£). 

Lemma II.7. Let r £ (0,1] and (3 £ (0,2] such that /3 — 1 + ?’ 2 / 0. We consider the 
equation 


1. Ifr 

2. Ifr 


cot (n6) + 


(2 - f3)r 


1 + r !: - !? eot ( 0)= 0 . 


r 2 + (/3 — 1 ) sin($) r 2 + (/3 — 1 ) 

1 , then ( 12 ) has exactly 2 n distinct solutions in [— tt, tt). 
|/3 — l| ^ { 0 , 2 }, we distinguish two subcases: 


( 12 ) 


(a) if f3 < 1, then (12) has exactly 2 n — 2 distinct solutions in [— 7 r, 7 r) for any 
n> 1, 

(b) if (3 > 1, then (12) has exactly 2n — 2 distinct solutions in [—7 r, 7r) if and only 
if n > 2 ( 2 -/ 3 ) • If n — 9 ( 2 —^ 3 ) > t/ien ( 12 ) has exactly 2 n distinct solutions in 
[—TT, 7t) . 

3. If j3 > 1, and ?’^{1, /3 — 1}, we distinguish two subcases: 


(a) if\<f3<l + r, then (12) has exactly 2 n — 2 distinct solutions in [—tt, tt) 

if and only if n > (1 ' % n - (i-r)(r+/3-i) > then ( 12 ) has exactl V 2n 

distinct solutions in [—tt, tt). 

(b) If (3 > 1+r, andn > 2, then (12) has exactly 2n —4 distinct solutions in [—tt, tt) 

if and only if n > ■ V (i+^-T-r-) < n ^ (i- r ^ r + +t3-i) ’ then ( 12 ) 

has exactly 2n — 2 distinct solutions in [—tt, tt). Finally, if n < 

f/ien (12) has exactly 2 n distinct solutions in [—tt, tt). 



4- If f3 < 1, and r ^ { 1, 1 — f3}, we distinguish two subcases: 

(a) if 1 — r 2 < f3 <\, then (12) has exactly 2n — 2 distinct solutions in [—7r, 7r) if 

and only if n> , 

(b) if 0 < /? < 1 — r 2 , we distinguish two subcases: 

i. ifO < fd < 1 — r, then (12) has exactly 2n — 4 distinct solutions in [— it, it) 

if and only if n > ■ If n < , then (12) has exactly 

2n distinct solutions in [—7T, n). 

ii. If 1 — r < (d < 1 — r 2 , then (12) has exactly 2n — 2 distinct solutions on 
[—7r, 7r) for any n > 2. 

5. If fd = 1 then (12) has exactly 2n — 2 distinct solutions in [—7r, 7r) if and only if 
n > jf-p. If n < jfp;, then (12) has exactly 2n distinct solutions in [—7r, it). 


III. MODEL SPACES AND MODEL OPERATORS 


This section lays down the theoretical framework on which our results are footed. We 
refer to [NN3] for an in-depth study of the topic. To the minimal polynomial m = niT of 
T G M n we associate a Blaschke product 


B{z) := n 

i 


Z - Xj 
1 - A iz 


The product is taken over all i such that the corresponding linear factor z — A* occurs 
in the minimal polynomial m. Thus, the numerator of B as defined here is exactly the 
associated minimal polynomial. The space of holomorphic functions on B is denoted by 
Hol{ B). The Hardy spaces considered here are 

1 r 2lx 

H -2 := {/ € Hol(B) : \\f\\ 2 H := sup — / |/(re^)| 2 d0 < oo}, 

0<r<l 27T J o 

and 

Hoo ■= {/ e Hol{ B) : ||/|| H := sup \f(z)\ < oo}. 

zeo 


Let B be the Blaschke product associated to the minimal polynomial of T € A4 n . We 
define the |m|-dimensional model space 

K b := H 2 © BH-2 := H 2 D {BH 2 ) L , 

where we employ the usual scalar product from L 2 {d B). The model operator Mb acts on 
Kb as 


Mb ■ I<b -a Kb 

f ^ M B (f) := P B {zf) 
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where Pb denotes the orthogonal projection on Kb- In other words, Mb is the compression 
of the multiplication operation by 2 to the model space Kb- Multiplication by z has 
operator norm 1 so that || Mb || < 1. Moreover, it is not hard to show that the eigenvalues 
of Mb are exactly the zeros of the corresponding Blaschke product, see [NN2]. 

For our subsequent discussion we will assume that T can be diagonalized. This does 
not result in any difficulties since the upper bounds obtained in the special case extend by 
continuity to bounds for non-diagonalizable T, see [OS]. One natural orthonormal basis 
for I\b (f° r diagonalizable T) is the Malmquist-Walsh basis {ek}k=i \ m \ with ([NN3, p. 
137]) 


with 


efc(z) 


(1 - |V| 2 ) 1/2 TT Z-A, 
1 1 A iZ 

i=i 


e i (z) 


d-|A 1 | 2 ) 1/2 

1 — Ai z 


IV. PROOFS 

A. Proofs of Theorem II.1 and Proposition II.2 

We recall [OS, Theorem III.2] that given a finite Blaschke product B = n /=1 i~p * z > °f 
degree N > 1, with zeros fi £ D, the components of the resolvent of the model operator 
Mb at any point (sC - {iq,..., I'/v} with respect to the Malmquist-Walsh basis are given 

by 


((C-Mb)-‘) 


l<i 5< 7<Ar 


To 


0 

1 n 

1 -ViC 11 s^l l-u a ( 

d-H 2 ) 1 ' 2 (1-F,I 2 ) 1/2 n t-^B 

i-viQ i-zqC •••*} i-^sC 


if i < j 
if i = j 
if i > j 


(13) 


We begin by proving Proposition II.2. 


Proof of Proposition II.2. We consider T € C n and assume for a moment that cr(T) C D. 
It follows from [OS, Theorem III.2] that for any ( £ C —cr(T) the resolvent of T is bounded 
by 


ll*(C, T)|| < \\R((, M b )\ , 

where B is the Blaschke product associated with iiit■ We conclude from (13) that 


((c- Mb ) 


r 0 


-1 


l<t, j<\m\ 


= < 


1 -AiC 1 -AjC 


if i<j 
if i = j 

nUfrff io>i 
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Taking the limit |A*| —>• 1 we have that for ( £ C — cr(T) the norm of the resolvent is 
bounded by the norm of the diagonal matrix 

/A 

l 

C-An 



which is dfTTJ- 


□ 


Proof of Theorem II. 1. We begin by showing that (5) holds. As before we suppose u{T) C 
D. The more general case cr(T) C B follows immediately by continuously moving the 
eigenvalue towards <90. As before we have that \\R{(, T)|| < ||i?(C, Mb) || and that 


((C - Mb) -1 ) 


l<i, jXlml 


r 0 if i < j 


(1-IA^) 1 / 2 (1 —jAjj^) 1 / 2 1 ■ . 

I i-AiC i-AiC ni- =J kx fe J 


where we abbreviated b\. = j- ^ . Recall that for any n x n matrices A = (ciij) and 
A = (ah), the condition |a^| < ah implies that ||A|| < ||A , ||. Hence, we estimate 


(c 


< < 


0 

j_ 1 

ii—Aid Ki 

(1-|6 A .| 2 )?(1-|6 A .| 2 )2 l 

m =j I bx k | 


if i < j 

if i = j 
if i > j 


where we made use of the well-known [JG] fact that 


j - |A| 2 1 I^a(C)| 2 

1-ACI 2 1-ICI 2 


By assumption min* |6^. | = p(f, a) = r and for any l with | b\ \ = r we have 


and since 


(C - Mb) 


< < 


i 1 —Aid r 0 

i-IMOl 2 


i-icr 


r i—j +1 


if i < j 

if i=j 
if i> j 


i-IMOl 2 = i-N 2 < 1 1 ~ N 2 < P 

1 - ICI 2 |1 - A;C| 2 “ d( l,dC) |1 - A;C| “ d(l,aCY 


with /3 = s(C, cr) we conclude that 



11 


(C - Mb);] 


< 


1 1 

d( 1, aQ rM 


0 

Jm|-1 

f3 r \m\-(i-j+1) 


if i < j 
if i = j 
if i > j 


Next we discuss the sharpness of (5). We consider ( £ [—1,1] and the model operator M B 
with B = bf and A £ (—1, 1). In this case d(l,a() = |1 — A£| and p((, a) = |&a(C)I and 
s((, a) = jj J^| . We find that 


0 if i < j 

M b )T 1 = - - * (&A(C)) n_1 if i = j . 

l ’ 3 aC) B ^ | ^^(&A(C)) n_(i_i+1) if i>j 

If b\(Q > 0 then the equality 

is clear. If h\(C) < 0 equality still holds because the Toeplitz matrices a = (ctfc)fc >0 and 
b = ((—l) fc dfc)fc>o have the same norm. To see this write D = diag(—l ,..., —1) and note 
that DT a D f = Tft. As D is unitary ||T a || = ||Tj,||. □ 


B. Proofs of Theorem II.4 and Propositions II.5, II.6 
Lemma IV. 1. For any ( £ ID) and r £ ( 0,1) 

7Z n (C, r) = n n (\(\,r). 

Proof. Indeed, writing f = |£|exp(i$) we have 11 — AC| = |l — exp(i#)A |d| and 

\\R((, T)|| = | (|C| — exp(—i0)T) _1 |, which entails that 

d( i, W) o me, t)i = d( i, ofr') ici) iNd ,toi 

with T' = exp (—i6)T. Clearly T £ C n O T' £ C n and p( cr(T)) = p(|C|, <r{T')). The 
same argument implies T £ T' £ Tff □ 

The lemma says that for our optimization task it is sufficient to consider real (. 

Proof of Theorem II.f. We consider a matrix T £ C n with p(f, er(T)) = r and the function 
s in (5). The function A i —> is maximized on the pseudo-hyperbolic circle of center 

ICI and radius r at A max = and takes the value /3 max = = i'd-^| ■ Thus 

using Lemma IV.1 and Theorem II.1 we find that TZ n {(, r) = IZ n (\(\, r) < ^ ||V r)( g max ||. 

According to the proof of Theorem II. 1 the model operator corresponding to Blaschke 
product B = bf with A £ (—1,1), A ^ |d satisfies 
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We choose A max = J^L _1 and have p(\(\, a) = r and s(|£|, a) = /3 max = 7 = 7^1 • It follows 
that lZ n {(,', r) > X lXj.^11 and the “worst” matrix for the extremal problem in (7) is the 
analytic Toeplitz matrix Mb with B = bf 

The monotonicity properties and the asymptotic behaviour of lZ n (C, r ) are consequences 
of our study of ||X r £ max ||, which is summarized in Proposition II. 6 . □ 

Proof of Proposition II. 5. We first prove (9). We consider T £ C n with spectrum a and 
assume for a moment that <72 C ID. It follows from [OS, Theorem III.2] that for any 
£ £ C — g the resolvent of T is bounded by 

me, t)| <|i?(c,M S )i, 

where B is the Blaschke product associated with P a . We conclude from (13) that the 
matrix of the resolvent of the model operator Mb at any point £ £ C — a (with respect to 
the Malmquist-Walsh basis) is given by a block-diagonal matrix 

Ai 0 \ 

0 A 2 )' 

where A\ £ _A/f|p CTi |(C) and A 2 £ Xf|p CT2 |(C) are respectly defined by (taking now the limit 
as \fi\ —>■ 1 for fi £ < 72 ): 


(A 


^'i<hj<\Pcr 


0 

1 

C-A i 

(1-IA.f) 1 / 2 (1-IAjj) 1 / 2 1—AfcC 

„ 1-AfC 1-AjC C-Afc 


if i < j 

if i = j 
if i > j 


and 


A‘2 


( 


1 

C—Mi 


V 


\ 

1 

C“Mn 2 ) 


where are the unimodular eigenvalues of Mb- To conclude, it remains to remark that 
the spectral norm of R {£, T) is nothing but 


max(|Ail , |A 2 |), 

that \\A -2 [j = dphrl)’ w hile fi 16 inequality 


\\Ai\\ < 


d (£, cn) 


X 


1, s(C,CTl) || ) 


follows directly from Theorem II.1. This proves the first inequality in (9). The second one 
is due to the fact that s(£,cn) is bounded from above by 2 and the equality refers to the 
resolvent estimate in [DS] or more generally to our Proposition II.6. 
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Now, let us prove (10). We consider T G C n with spectrum a = G\ U a 2 , n\ = \P ai | and 
ri 2 = \P a2 1 such that max^d |A| < p\. We fix f € c© — a. The last inequality, combined 
with the trivial observation that s(£,ci) < pi, proves that 

To complete the proof of (10), we first remark that by rotation invariance (as in the proof of 
Lemma IV.1) the sup defined in the statement of (10) does not depend on ( and coincides 
with its value at point ( = 1. We consider the spectrum a = G\ U cr2 defined by setting 
Pen = (z — p i) ni and P a2 = (z + l)™ 2 . Let B be the Blaschke product associated with P a . 
We repeat the proof of (9) to the model operator Mb- The matrix of its resolvent at point 
1 (with respect to the Malmquist-Walsh basis) is given by the block-diagonal matrix 

Aj 0 \ 

0 A 2 J’ 

where 


1 


1 


(Ai)i<ij< m - d{h ai)B{X) 

_ ^1,1+pi 


0 

(Wi)) 


n—1 


if i< 3 
if i = j 


L ^(^(l)) n - ( ^' +1) if i>j 


I 1 - Pi I 

and A -2 = where /„ 2 is the identity matrix of A4 n2 (C). Thus, 


||-R(1, M b ) I = max( 


\l-pi\ ^ Xl ’ 1+P1 ^ ’ 2^’ 




but d( 1, a) = d( 1, (J\) = |1 — pi|, and as a consequence: 

d{ 1, cr) ||i?(l, M b ) I = max(||Xi 1 i+pi I 

The above maximum is equal to ||Xi i+ Pl || : indeed, according to Proposition II.6, 

1 


II^M+pJ > ll^i.ill = 


2 sin( 


4ni+2 > 


and the latter quantity 


1 


2 sin( 


4n 1 +2 


is always greater than i (increasing with n\) since 


m > 1. The last part of the statement relies again on Proposition II.6 and the fact that 
P^\\X 1 , i_l_p | is increasing on [0, 1]. □ 


Proof of Proposition II. 6. Our approach is guided by the methods developed in [EE, DS, 
OS]. Often it will be convenient to transform the Toeplitz matrix X r ^p to an associated 
Hankel matrix X r ^p, see [DS, OS]. We set X r ^p = X r ^J, where 


J = 


/ 

\1 
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X is a Hermitian Hankel matrix that satisfies ||Af|| = \\X r «||. Hence, in our discussion we 
can focus on the largest in magnitude eigenvalue of X. (We omit indices in X for simpler 
notation.) 

%) In the case f3 = 1 — r 2 we have 


1 _1 
A y R — 
y>n ' ’ r y>n 


0 if i < j 

rn_1 . . Xi=3 =( M b« r ) l< M <n. 

(1 — r 2 )r n J+1 ) if i > j 


where the second equality is (13). It is well known [NN1, Theorem 3.12, p.147], [OS, 
Lemma III.5] that (see the cited articles for a discussion) 


Mb 1 


I^/zIHoo/BHoo 


1 


n 

AEcr 


i 

lA|' 


In particular, we find 



frr and \\X r ,p\\ = 1. 


ii) X is an entry-wise non-negative matrix. A standard result in Perron-Froebenius 
theory [HM, Theorem 4.2] asserts that the spectral radius of such matrices is a non-negative 
eigenvalue and the corresponding eigenvector is entry-wise non-negative 1 . Furthermore the 
spectral radius is monotonically increasing in the matrix-entries [HM, Corollary 2.1], i.e. 
X < B entry-wise implies p(X) < p(B). 

in) Here, we show a simple proof for f3 > 1. The more general discussion under the 
condition 1 — r 2 < j3 can be found below. Let X n G A4 n and X n+ i G M n + 1 - We show 
||X n || < |A n+ i||. Let x* with ||x*|| = 1 and X n x* = ||X n ||a;*. Computing entries of vectors 
and using (3 > 1 we see that 


X„x* 


< 


X 


n +1 


iv) We turn to the case r 2 +/3 —1 ^ 0. The discussion is guided by [EE]. The eigenvalues 
of X 2 are the eigenvalues of X squared. Hence, we are looking for the largest A 2 such that 

0 = det(X 2 - A 2 ) = det(X - A) det(X + A). 

The matrix X 2 — A 2 can be computed explicitly. We display this matrix for n = 5 for 
simplicity, see Appendix V for the general matrix. If n = 5 we have 


X 2 - A 2 = 


f A 2 +/3 2 ±^+rS 

^+Pr 7 

P 2 r 2 ^-+pr 6 

P 2 r 34 -f^-+Pr 5 

Pr 4 \ 

P 2 r^+pr 7 

A 2 +/3 2 r 2 lp^ + C 

P 2 r 3 ^-+Pr 7 

P 2 r 4 ±=^+pr* 

Pr 5 

P 2 r 2 ^-+Pr 6 

P 2 r 3 ^+pr 7 

A 2 +A 4 ^+r 8 

P 2 r 5 ^+Pr 7 

pr 6 

P 2 r*±f^+pr 5 

P 2 r 4 ±=^+pr e 

P 2 r 5 ^+Pr 7 

\2 + p2 r 6l^_ +r S 

Pr 7 

\ Pr 4 

fir 5 

1 3r e 

pr 7 

—A 2 + r 8 / 


The determinant of the above matrix is transformed by the following steps. 


1 In fact, as can be seen easily by computing powers, X is primitive. 
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1. We divide the k th column and the k th row of ( X 2 — A 2 ) by r k 1 . 

2. In the determinant resulting from (1), we subtract the k th column from the (k— l) th 
and leave the n th unchanged. 

3. In the determinant resulting from (2), we subtract the (k — \) th row from the k th 
and leave the n th unchanged. We obtain the determinant of a tri-diagonal n X n 
matrix. 


4. In the determinant resulting from (3), we multiply the k th row by r 2k 1 
We obtain that det(W 2 — A 2 ) = 0 iff 



/ 7 

0! 

0 

••• 0 \ 


( 7 

k 

0 

••• 0\ 


a 

7 

0! 




7 

k 


0 = det 

0 

a 


0 

= det 

0 

k 


0 





7 a' 





7 k 


Vo 


0 

ay) 


u 


0 

k y / 


where a = A 2 r + r 2n 1 (/3 — 1), a' = ^,y = — r\ 2 + r 2n 1 and 7 = — A 2 (r + 1/r) + r 2n 2 (r + 
i(/3 — l) 2 ) and we set y = a/r in the last step. We have that 


/ 7 k 0 ■■■ 0 \ 


det 


f* 7 M : 
0 n ■ • ■ • 0 


■■ 7 M 

V 0 • • • 0 Li y J 


k 

0 ••• 0 \ 


( 7 k 0 • • • 

0 ^ 


7 

k : 


k 7 k 



k 

0 

+ det 

0 y 

0 



■ • • 7 n 


: 7 

0 



0 k 7 / 


O 

O 

y-7 ) 



( 7 k 0 

••• 0 \ 


( 7 k 

0 ••• 0 \ 


k 7 k 



k 7 

k : 

det 

0 y 

0 

+ {y- 7 ) det 

0 y 

0 



7 k 



■ • • 7 y 


V0 ••• 0 

k 7 / 


U ••• 

0 y 7/ 


where in the last sum the first determinant is of order n, while the second is of order n — 1 . 
The assumption /x = 0 is equivalent to A 2 = r 2n ~ 2 ( 1 — (3) and the determinant equation 
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turns into 7 y = 0. In this case we obtain a solution if (3 = 0 or r 2 + f3 — 1 = 0 but these 
choices of (3 are treated differently. (Note that A = 1 implies [i 7 ^ 0.) For the determinant 
of a tri-diagonal Toeplitz matrix we have the well known formula [MV, EP, EE, OS] 


a 1 0 

1 a 1 


■. 0 — U n (a/ 2 ), 


: , - . 1 
V 0 • • • 0 1 a) 

where U n denotes the n-th degree Chebychev polynomial of the second kind. Hence the 
eigenvalues of X correspond to the roots of the polynomial equation 

V n+ 1 u n (+ i-i n (y - 1 ) U n —1 = 0 . 


v) We make use of properties of U n to reveal properties of the roots of equation (11). 
Let us consider A 2 so that there is real 9 with 7/2 fi = cos (9), i.e. 

\2 = 2(n—1) W ~ l ? + r2 ~ 2r (P ~ 1 ) cos (#) (u) 

1 + r 2 + 2 r cos (9) 

Using U n (cos(9)) = Equation (11) turns into 

sin((ra + 1)9) 1 A 2 — r 2n ~ 2 ((3 — l ) 2 sin(n$) 

sin($) r A 2 + r 2n ~ 2 (/3 — 1) sin(0) ’ 

which plugging back A 2 is 

Similarly, if r j/2/a G [l,oo) we set 'y/2/a = cosh($) with 9 G [0, 00 ), which leads to [MH, 
Section 1.4] 

\ 2 = 2 (n— 1 ) (/? - l ) 2 + r 2 - 2r(f3 - 1) cosh(fl) 

1 + r 2 + 2r cosh($) 


sinh ((n + 1)9) _ (2 — /3)r + 2(1 — /?) cosh(0) 

sinh (■ n9 ) r 2 + (3 — 1 

Finally, in case that 7 / 2 /t G [—1, — 00 ) we set 7 / 2 [i = — cosh(0) with 9 G [0, 00 ) and recall 
that U n (—x) = (—1 ) n U n (x) [MH, Section 1.4] to find 

A 2 = r 2(n-i) (ff ~ !) 2 + r ' 2 + 2r (fi ~ x ) cosh( 6 >) 

1 + r 2 — 2 r cosh($) 


( 18 ) 
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with 


sinh ((n + 1)9) (2 — f})r — 2(1 — f3) cosh(#) 


sinh ( nd ) 


r 2 + /3 — 1 


(19) 


We consider the first case. Any solution 6/ ; of (15) leads to a solution A| of det(A' 2 — 
A 2 ) = 0 via (14). An important point about Equation (14) is that all such solutions || 2 
can be bounded by a function that decays exponentially with n. This allows us to study 
properties of ||X|| by counting solutions of equation (15). We observe that 


1 


sin(0) 

1 


-T Too 

as 9 —> 0 + 

and 

tan(0) ^ + °° 

-T Too 

as 9 —> tt~ 

and 

tan(0) ^ °° 


as 




sin(0) 

It follows that ifl — r 2 < (3 < 1 + r then (see Lemma II.7 for details) 

1 


as 


1 r 2 - (/3 — 1) (2 - (3)r 


tan(0) r 2 + (3 — 1 
1 r 2 — (/3 — 1) 


+ 


r 2 + (3 — 1 sin(0) 

(2 - P)r 1 


+oo as 9 —> 0 + 


tan(0) r 2 + (3 — 1 r 2 + (3 — 1 sin(0) 

As — cot (nd) has n branches in the interval [0,7r] and since 

— cot(n0) -T —oo as 6 —> 0“ 

— cot (nd) —> Too as 


Too as 


7r 


7T 


7T 


we conclude that Equation (15) has n — 1 solutions 9j. G (0, tt) (if n is large enough, see 
Lemma II.7). The corresponding eigenvalues A& of X approach 0 exponentially fast, so 
that for the largest eigenvalue of X we have A* Tr(X). On the other hand direct 
consideration of the trace shows that at fixed (3 and r we have 

TV(jn -r —as n — > oo. 

1 — v 2 

If f3 > lTrwe have only n — 2 real solutions 9j~ to Equation (15), see Lemma II.7. We 
consider Equations (16), (17) and Equations (18), (19). As n gets large quickly 

converges to e e with 9 G [0, oo). By the intermediate value theorem (17) has a solution for 
f3 > IT r and sufficiently large n. Any A 2 as in Equation (16) decays exponentially with 
n. Hence, in case /3 > 1 T r we obtain n — 1 eigenvalues that decay exponentially with n, 
as before. In the asymptotic limit of large n we find that = \/r becomes a solution of 
(19) and A* in (18) does not decay exponentially with n, which identifies the latter as the 
largest in magnitude solution for large enough n. Hence, we have that ||X|| -T f3/( 1 — r 2 ) 
as n —> oo. 

in) (second part) We have already seen that || X r jj || grows monotonically with n if /3 > 1. 
Hence, we consider 1 — r 2 < f3 < 1. As outlined in v) there are for large enough 2 n, see 


2 When there are n solutions to (15), monotonicity of |X r)J g| can be readily be seen in (14) and (15). 
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Lemma II.7, n— 1 solutions 9 *. to Equation (15) and the corresponding eigenvalues, see (14), 
can be bounded as 

r 2 '”- 1 ’ (l ~ YV7) 2 - M " T^) 2 • (2D 

We write X^ and X^ n ~^ for the matrix X in n and n — 1 dimensions. Let || -||^ denote 
the Frobenius norm. Direct computation shows that 

X(n) 2 _ x(n- 1 ) 2 = nr 2 (n-l) + _ 1 )^ r 2 (n- 2 ) _ , _ fo r 2 (n- 2 )_ 

FR FR 

On the other hand it is well known that \ 2 FR = YPi= 1 4 2 , where the A -”" 1 are the 

eigenvalues of X^- n \ Writing A*’( n ) for the largest in magnitude eigenvalue of X^ and 
making use the mentioned properties and the estimates ( 20 ) we find that 

|^*,(n )|2 _ 1)|2 

> nr 2(n_1) + (n - 1 )/3 2 ?’ 2(n-2) - (n - l)r 2(n_2) 

- (n- l)r 2( ” _1) (l - y~^) 2 
+ (n-2)r 2 ^ 

\ 1 + rj 

The lower bound is a degree-2 polynomial in /3, which is convex if r € (0,1/2). Hence, 
the lower bound is an increasing function of /3 € [1 — r 2 ,1]. As /3 = 1 — r 2 implies 
|^*,(n )|2 _ |^*,(n-i )|2 _ q proves monotonicity in n if r € (0,1/2). The range for r can 
be extended by using better bounds in (20). Here we note that at the endpoints of the 
interval for f3 we already have established monotonicity for any r. Furthermore the norm 
is a monotonous, convex functions in /?. Any of the branches of — cot (n9) is convex and 
the terms 

AAAiL and ( 2 ^ 

r 2 + f3 — 1 r 2 + /3 — 1 

in (15) are decreasing in € [1 — r 2 ,1]. It follows that the solutions 8^ of (15) satisfy 

O^iP 0 ) - e^(Pi) < ^ n-1) (/3o) - ^ n_1) (/3i) 

for Po < Pi- In conclusion the quantity 

nr 2 (n-l) + _ l )/ 3 2 r 2 (n- 2 ) _ („ _ 1 )r . 2 (n- 2 ) + £ | A £—D |2 _ £ | A jA )|2 

k =1 k =1 

is increasing in /3. 

vi) is a direct consequence of in) and v). □ 
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Proof of Lemma II. 7. For fixed (3 G [0, 2], r € (0, 1] we for now focus our attention on real 
values of 8 and write shortly 




r 2 -((3- 1) 
r 2 + (/3 — 1 ) 


cot( 0 ) + 


(2 - p)r 1 
/3 — 1 + r 2 sin(0) 


such that Equation (12) becomes 


l Pr,/3(8) = —cot(nd), 8 G [— 7 T, it). 


We start with some elementary observations. 


1. We restrict to the interval [0, ir) as all functions are symmetric with respect to the 
point ( 7 r, 0 ). 

2 . cot(n 0 ) has poles in [ 0 , 7 r) at 8 = ^ 7 r, k = 1 , ...,n — 1 . 

3. The functions sin ^ and cot(0) have poles at 0 and 7 r, cot(0) ~ as 6 tends to 
0 + and cot( 0 ) ~ — sin * ^ as 8 tends to 7 r~. 

4. Looking at the polynomials: 

P 0 (r) = r 2 - (2 - (3)r - (/3 - 1 ) = (r - 1 )(r - (1 - /3)), 

and 


Qp{ r ) = r 2 + (2 - /3)r - (/3 - 1 ) = (r + l)(r - (/? - 1 )), 


we have that Pp{r) < 0 for r € (1 — /?, 1) and Pp(r) > 0 elsewhere, while Qp(r) < 0 
for r G (—1, /3 — 1) and Qp(r) > 0 elsewhere. 

5. Taking (c) and (d) into account, we have 




r 2 + (2 — /3)r — (/? — 1 ) 1 

(3 — 1 + ?’ 2 sin(0) 


<3/?(r) 1 

(3 — 1 + ?’ 2 sin(0) ’ 


while 


Pr,p{0) 


r 2 - (2 - /3)r — (/3 — 1) 1 

/3 — 1 + ?’ 2 sin(0) 


-Ppir) 1 
/3 — 1 + ?’ 2 8111(0) ' 


We now prove the points in the lemma. 

1. The case r = 1. In this case, the structure of (12) is discussed in [OS, Section B], 
We now suppose that r € (0, 1). 


2. The case r = \f3 — 1| ^ {0, 2}. 
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(a) The subcase j3 < 1. We have that 


<Pr,p(&) 


(—2 + P) (cos(0) + 1) 
sin( 0 ) p 


and as a consequence, p r ,/3 < 0 on ( 0 , 7 r), limg_ > ,o+ p r ,p ( 0 ) = — oo, while 
liiri fl _ {p r j g(0) = 0. Thus, the graph of p r ,p intersects n — 1 times with 
the graph of 9 i->- — cot(n0) on [^, 7 r). Moreover, 

p r ,p{6) (2 — P) sin(0) sin(n0) 

— cot(n 0 ) /3cos(ra0)(l — cos( 0 )) ’ 

and in particular, for any 0 € ( 0 , ^), ^cot(nd) > which can be seen by 

computing derivatives with respect to 9: 9 H > _Tcot(ne) is increasing on (0, ^), 
with 


<Pr,p{0) _ 2n(2 — P) 


lim 

0 —>o+ —cot(ntf) 


P 


> 2n > 1 . 


The graph of p r ,p does not intersect the one of 0 h-> — cot(n0) on [0, #-). We 
conclude that in this case (12) has 2n — 2 real solutions 9k in [— 7 r, 7 r). 

(b) The subcase P > 1 . We have that 


Pr,p{Q) 


(-2 + P) (cos (9) - 1) 
sin( 0 ) p 


and as a consequence, p r ,p is increasing on (0, 7 r), p r ,p > 0 on (0, 7 r), 


lim 


Pr, /?(0) — 0j 


while lini fl _ ip r p(9) = +oo. Hence the graph of p r ,p intersects n — 1 times 
the one of — cot(n0) on (0, ]. Moreover, 

ip rj p(9) (2 — P) sin(9) sin(n9) 

— cot (n9) /3cos(ra0)(l + cos(0)) ’ 

and in particular, for any 9 € ( ^ 2 ^ > 7r )> -cot(n0) - > 2n ^ 13 " 1 using the fact 
(as above) that 9 1 —>• ^cot/nO) ti me decreasing on (^— 5 ^—, 7 r), with 

li m _ 2n(2 - P) 

0 -» 7 T— — cot(n 0 ) P 

Adding the positivity and the continuity of 9 i-> ^cot/nfl) 011 ^~ 2 rP~' 7r )’ aR d 


lim 

6 — > (^ n ~i) 7r 


— cot(n 0 ) 


= +00, 
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the intermediate value theorem guarantees that the graph of <p r g does not 
intersect the one of 6 i->- — cot (nO) on ( ^2 ~n ^ > 7r ) if and on ly if 2 "^j ^ > 1 
(i.e n > 2 (<f- ft) )• We conclude that in this case (12) has 2n — 2 real solutions 
9k in [—7T, 7r) if and only if n > 2 ( 2 -fi) • 

From now on, we suppose that r 7^ 1 — (3, and r 7^ f3 — 1. Let us first notice that in 
case /? — 1 + r 2 7^ 0, we have 


li m ^(A) = n(l + r)((3 - 1 - r) 
0^0+—cot (n0) /? — 1 + r 2 


and 


li m = w(l -r)(r + / 3 - 1 ) 

0-771— — cot(n 0 ) /? — 1 + r 2 

Below, we repeat the reasonings of (2) for the remaining cases (3), (4) and (5). 

3. The case (5 > 1, and r ^ {1, /? — 1}. With (e), Pp(r) < 0, for any r £ [0, 1). In 
particular, for any r € [0, 1), 


lim ( p r ,p{0 ) = +00, 

0 — >7T 


while 


lim (p r ,p(0) 
0—>O+ 


—00 if r £ [0, j3 — 1) 
+00 if r £ (f3-l, 1) 


because of the property of Qp(r) in (d). 


(a) The subcase 1 < j3 < 1 + r. In this case, ip r , p > 0 on (0, 7r) and 

lim ip r a(9) = lim <p r p(9) = +00. 

0—70+ e->7r- 

Moreover, ip r « has a global minimum at some point 9q £ (0, 7r) such that 
(p r p is decreasing on (0, $o] and increasing on [9q, it). Thus, the curve of 
intersects the one of 0 1—>■ — cot(n@) exactly n — 1 times on (0, ( 2n ~ 1 ) 7r ] i 
On the interval 7r), the curve of <p r ,p does not intersect the one of 

9 1—^ — cot(nfl) if and only if lim^^- _fcot(nfl) > b (else, the number of real 
intersections on [0, it) is n). Indeed as before, 9 1—>• is (positive, con¬ 
tinuous and) decreasing on 7r). We conclude that in this case (12) 

has 2n — 2 real solutions 9k in [— tt, 7r) if and only if n > j • 
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(b) The subcase /? > 1 + r. In this case, tp r ,p is increasing on (0, tt), and 
lim 0 _ ) , o + (p r ,p(6) = — oo, while lim ip ri p(9) = +oo. Thus, the curve of 

tp r ,p intersects the one of 6 i —> — cot (n9) exactly n — 2 times on [T-, 
Moreover, on (0, U 7 r), 6 i —y ^^ot(n9) Positive and continuous; it 

is increasing on (0, -ff~) and decreasing on p—, tt). Thus, the curve of <p r ,f) 
does not intersect the one of 6 i —> — cot (nd) on ( 0 , ff-) (resp. on tt)) 


if and only if lim^ 0 + > 1 (resp. lim^- T C ot(ng) > x )- We conclude 

that in this case (12) has 2n — 4 real solutions 6k in [— 7 r, tt) if and only if 

0—1+r 2 
(1—r)(r+/3—1) ’ 


“Pr, /?(^) 


/3-1+r 2 4~l+r 2 , 

n > (1—r)(r+/3—1) ’ smce (l+r)(/3—1—r) < 


4. The case j3 < 1 and r ^ {1, 1 — /?}. With (e), Qp(r) > 0, for any r G [0, 1). In 
particular, for any r G [0, 1), 


lim <p r ,p(6) 

0->O+ 


+oo if /? — 1 T r 2 > 0 
—oo if /? — 1 T r 2 < 0 


while 


lim <^( 0 ) 

0— >-7T 


Too if /3 — 1 T r 2 > 0 

Too if [3 — 1 T r 2 < 0 and r € [0, 1 — fi) , 

—oo if j3 — 1 T r 2 < 0 and r G (1 — /?, 1) 


because of the property of Pp{r) in (d). 


(a) The subcase 1 — r 2 < /3 < 1. In this case, <p rt p > 0 on (0, i r), and 

limg_ >0 + <p r ,p{6) = lim ip r ^(6) = Too. Moreover, has a global min¬ 
imum at some point 0q G (0, n) such that (p r ,p is decreasing on (0, do] and 
increasing on [6o, n). We conclude as in (3)-(a) that in this case, (12) has 
2 n — 2 real solutions 9k in [— tt, n) if and only if n > . 

(b) The subcase 0 < (5 < 1 — r 2 . 


i. If 0 < j3 < 1 — r, (p rt p is increasing on (0, tt), lim fl + <p r ,p{ 6 ) = —oo, and 

linig^j.^- y>r,p{ 6 ) = Too. We conclude by a similar reasoning as the one in 

( 3 )-(b) that in this case also ( 12 ) has 2 n — 4 real solutions 6 k in [—tt, tt) if 

and onlv if n ~> / 3 ~ 1+ ' r2 s ; 1lrP p-l+r 2 , p-l+r 2 

and only it n> ( 1 _ r .)( r+/3 _ 1 ), since ( 1+r )( /3 _ 1 _ r ) < ( 1 _ r )( r+( a_ 1 ) ■ 


ii. If l — r<(3<\ — r 2 , <p r ,p(6) < 0, \/6 € (0, tt), and lim 6 i_ i>0 + (p rj p(9) = 
lim fl <p r ,p(9) = —oo. Moreover, has a global maximum 9q G (0, tt) 
such that <p r ,p is increasing on (0, 0q] and decreasing on [0q, tt). Thus, the 
curve of tp r ,p intersects the one of 6 H >■ — cot (nO) exactly n — 1 times on 
[■^~, tt). in the interval (0, ff-), the curve of <p r ,p does not intersect the one 

of 6 i— y — cot (n0) if and only if lim^ 0 + ^he latter condition 

is always satisfied since n ^ 1+ ^f +r l — = > n. We conclude 

that in this case, ( 12 ) has 2 n — 2 real solutions 6k on [—tt, tt). 



5. The case /3 = 1. In this case, 
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<Pr,p(P) 


1 + ?’COs(0) 
rsin(0) 


and thus (p r ^(0) > 0, MO G (0, 7 r), and lining +ip r ^(0) = lim^ ^ <p ri p(9) = +oo. 
Moreover, (p r ,p has a global minimum at some point 6q G (0, 7 r) such that 
is decreasing on (0, $o] an d increasing on [0 q, 7 r). In particular, the graph of <^v i( g 
intersects the one of 6 H» — cot(nO) exactly n — 1 times on (0, On the 

interval (, 7r), the function 6 * i-G is decreasing and thus the curve of 

(p r ,/3 does not intersect the one of 6 1 —» — cot(n 0 ) if and only if lim^^- ^ot(n 0 ) > 1 - 
We conclude that in this case, (12) has 2 n — 2 real solutions 0 on [—7T, 7 t) if and 
only if n > (else (12) has 2n real solutions in [— 7 r, 7 t)). 


□ 


V. APPENDIX 

The columns C %, i = 1,..., n of the n x n matrix A — A 2 are: 


Oi 


— A 2 -{- /3 2 (1 -)- ? >2 -)-••• T 2) j _j_ j,2.(n—1) 

(3^ T (1 T ? >2 T T 2 ‘^ n ~) -f- f3 i ) — 1 
/ 3 2 r 2 (l + ?’ 2 + y,2.(n— 4 )^ _|_ p j,2.(n—1)—2 


f3 2 r n 3 (1 + r 2 ) +/3 r n+1 

f3 2 r n -2 + p r n 
P r n_1 


C 2 


(3 2 r( 1 T ? >2 T T 2 '^ n ~ 3 ^ -{- i ) — 1 

—A 2 T /? 2 r 2 (l -)- r 2 7-2.(n— 3 )^ _j_ ^.2.(n—1) 

P‘2 ^3^ _j_ ^,2 _j_ j,2.(n 4 )^ _|_ p pi.(n 1) 1 
/3 2 r n-2( 1 + r 2) +/3r n+2 

/3 2 r n_1 + /3r n+1 
(3 r n 



Cs 


(3^ r 2 (l + r 2 + ^ 2 .(n—4)^ _|_ y,2.(n— 1)—2 

^2 t»3^ -|- 7*2 _|_ y,2.(n 4)^ _|_ ^ j,2.(n 1 ) —1 

— A 2 + /3 2 r 4 (l + r 2 + . . . + 7 " 2 .(n— 4 )^ _|_ ^ 2 .(n— 1 ) 

^2 ^5^ -|- 7*2 _|_ ^ _|_ y,2.(n 5)^ _|_ ^ y,2.(n 1 ) —1 


/3 2 r n 1 (1 + r 2 ) +/3 r n+3 
(3 2 T n + (3 r n4 ” 2 
(3 r n+1 


^2 ^3 ^ _|_ ^2 _|_ j,2.(n 5) ^ _|_ ^ j,2.{n 1) 3 

/3 2 7^(1 + 7" 2 -(- 7*2.(n—5)^ _|_ 0 y,2.(n—l)—2 
^2 7*5^2 -(- 7*2 _|_ ^2.(n—5)^ _|_ ^ y,2.(n— 1) 1 
_^2 ^2 r 6^ _|_ 7*2 _|_ # _ _|_ 7,2.(n—5)^ _|_ r 2.(n-l) 

/^ 2 7^(1 + 7^ 2 7*2*(^ — 6)^ _|_ 0 j,2.(n—l) 1 

/3 2 r n (l + r 2 ) + /? r n+4 

0 2 ^n+1 _|_ ^ ^72+3 

/3r n+2 


(3 2 r n 4 (1 + ?’ 2 + ?’ 4 ) +/3 r n+2 

0 r n- 3 ^ _|_ ? .2 _|_ ? , 4 ^ _|_ p j.n+3 

f3 2 r n ~ 2 ( 1 + ?’ 2 + ?’ 4 ) + /3 r ,n+4 


02 02n 10^ _j_ j,2 _|_ ^,4^ _j_ j,2n 4 
02 r 2n—9 0 _|_ ^2 _|_ r 4^j _|_ r 2n—3 

—X 2 + /3 2 r 2n ~ 8 ( 1 + r 2 + r 4 ) + ? ,2n ~ 2 
/3 2 r 2n ~ 7 (l -|- r 2 ) -|- f3 r 2n-3 
0 ? ,2n-6 _|_ ^ ? ,2n-4 

/3 r 2n “ 5 


f3 2 v n 3 (1 + ?’ 2 ) +/3r n+1 

0 r n-20 _|_ ? ,2j _|_ p r n +2 


Cn—2 


0 r 2n-90 _|_ r 2j _|_ p r 2n-5 
0 r 2n-%0 _|_ ? ,2j _|_ p 0U-4 
0 r 2n-7 0 _|_ r 2j _|_ p r 2n—3 

—A 2 + /3 2 r 2n_6 (l + r 2 ) + ? ,2n ~ 2 

0 r 2n-5 p ? ,2n-3 
p r 2n-4 
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p 2 r n - 2 + pr n 

p2 r n—l p r n+l 

p 2 r n + f3r n+2 


and 


C n -1 = 


P 2 r 2n-5 _|_ p ? ,2n-3 
— \2 _|_ ? ,2n—4 _|_ r 2n—2 

/3r 2n ~ 3 


f3v n 1 
j3r n 
j3r n+l 


C n = 


f3 r 2n ~ 4 

P ? ,2n-3 

—A 2 + r 2n_2 
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